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The orbital dynamics of the binary point-mass systems is ambiguous at the third post-Newtonian 
order of approximation. The static ambiguity is known to be related to the difference between the 
Brill-Lindquist solution and the Misner-Lindquist solution of the time-symmetric conformally flat 
initial value problem for binary black holes. The kinetic ambiguity is noticed to violate in general the 
standard relation between the center-of-mass velocity and the total linear momentum as demanded 
f*"^ , by global Lorentz invariancc. 
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I. INTRODUCTION 

The motion of binary systems is one of the most important problems in general relativity, particularly in view 
of the future detection of gravitational waves from such systems. The simplest two-body problem is the one where 
the components of the system are objects which are as spherically symmetric and point-like as possible. In general 
^ ' relativity, those objects are black holes. In spite of the extension of black holes, in binary orbit approximate calculations 
of post- Newtonian (PN) type (which use expansions in powers of 1/c, where c denotes the velocity of light), the ansatz 
of Dirac delta functions to describe point-like sources turned out to be very successful up to the 2.5PN approximation 
and at the 3.5PN order of approximation too ||. At the 3PN order, i.e. at the order (1/c 2 ) 3 , only partial success 
was achievable Two terms in the binary point-mass calculations came out ambiguously: a kinetic one, depending 
on the bodies' linear momenta squared and a static one, not depending on the momenta j5j. The authors were 
able to show that the static ambiguity is intimately related to the fact that different exact solutions (Brill-Lindquist 
and Misner-Lindquist) of the time-symmetric and conformally flat initial value problem for two black holes do exist, 
which at the time reveal the emergence of black holes in binary point-mass calculations at the 3PN order |J. In 
this paper we confront the kinetic ambiguity with the center-of-mass motion and claim that it can be fixed through 
the standard relation between the center-of-mass velocity and the total linear momentum which results from global 
Lorentz invariance. 

We employ the following notation: x a = [x l a ) (a = 1,2; i = 1,2,3) denotes the position of the ath point mass in 
the 3-dimensional Euclidean space endowed with a standard Euclidean scalar product (denoted by a dot). We also 
define r i2 = xi — x 2 , r i2 = |r i2 1 , ni2 = i"i 2 /ri 2 ; | • | stands here for the Euclidean length of a vector. The linear 
momentum of the ath body is denoted by p a = (p a i). 
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II. THE GENERALIZED CONSERVATIVE HAMILTONIAN TO 3PN ORDER 



The post-Newtonian approximation order 2.5PN (as well 3.5PN) is of purely dissipative character. Subtraction 
of the 2.5PN terms from the 3PN Hamiltonian results in the conservative 3PN Hamiltonian. This Hamiltonian is 
primarily of higher order as it depends on time derivatives of the positions and momenta of the bodies Q . 

The conservative 3PN higher-order Hamiltonian for two-body point-mass systems was calculated in Ref. j|. We 
present here the explicit form of this Hamiltonian in the center-of-mass reference frame (where pi + p 2 =0). It is 
convenient to introduce the following reduced variables 

X i -X2 Pl = _P2 t= J_ £NR = £^ (1) 

~ GM ' n M ' GM' m ' W 

where 

mim 2 fi mim 2 . , 

M M (mi + m 2 ) 
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In Eq. ([j]) the superscript NR denotes a 'non-relativistic' Hamiltonian, i.e. the Hamiltonian without the rest-mass 
contribution Mc 2 . We also introduce the reduced-time derivatives: 



. _ dr d(xi — x 2 ) . _ dp G dpi G dp 2 
dt dt P dt v dt v dt 

The reduced two-point-mass conservative Hamiltonian up to the 3PN order reads 
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ff WK (r, p, r, p) = ff N (r, p) + -H 1PN (r, p) 



where (here r = |r| and n = r/r) 



1 - , , 1 - 
— i?2PN (r, p) + — iJ 3 PN (r, p, r, p) , 



~ / s P 2 1 

H N (r,p) = — -~, 



(3) 



(4) 



(5) 



£iPN (r, p) = i(3i/ - l)(p 2 ) 2 - 5 [(3 + ^)P 2 + i/(n ■ p) 2 ] J + 



(6) 



£ 2PN (r, p) = — (1 - 5v + hv 2 ) (p 2 ) 3 



-~ [(5 - 20^ - 3^ 2 ) (p 2 ) 2 - 2^ 2 (n ■ p) 2 p 2 - 3^ 2 (n ■ p) 4 ] I 
-i[(5 + 8,)p 2 + 3Kn- P ) 2 ]^-i(l + 3,)l 



(7) 



£ 3PN (r, p, r, p) = — (-5 + 35^ - 70^ 2 + 35^ 3 ) (p 2 ) 4 



,2\3 
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+— { (-7 + A2v - 53z/ 2 - 6i/ 3 ) (p 
16 I 

+ (1-2*V [2(n.p) 2 (p 2 ) 2 + 3(n.p) 4 p 2 ]}i 

3 (-27+ 140^ + 96j/ 2 ) (p 2 ) 2 
+6(8 + 25v)u(n ■ pfp 2 - (35 - 2Q7v)v{n ■ pf 
+ — — { [-4800 - 3(8944 - 315tt 2 )^ - 7808i/ 2 ] p 
+9 (2672 - 315tt 2 + 448j/) i/(n ■ P) 2 }^ 
12 + (872 - 63tt 2 )i/ 



1 r 



-Z?(r,p,r,p) + fi(r,p). 



(8) 



In Eq. (||) £> (r, p, f , p) denotes that part of the 3PN Hamiltonian i/3PN which depends on the time derivatives r 
and p. Its explicit form reads 



2 



1 12 



4(n-p)(n-p)(p-p)-5(n-p) 2 p 2 
5(n • p) 2 p 2 - (n • p) 2 (n • p) 2 + 13p 2 p 2 + 2(p • p) 2 



5 = ^ 2 {(n.r)(p.p)[5p 2 + (n.p) 2 ] 

-(n • p) [p 2 (p • r) + 2(p • p)(p • r)] - l(n ■ p) 3 (p • r) 
+ (n.p)[p 2 + (n.p) 2 ] [(n.p)(n.r)-(p.r)]}, 



-i = ^{ 2 ( 17 - 10 ^)( n ' P)( n ' P)( n ' r) - (15 + 22^)(n • p) 2 (n • 
-(51 - 8i/)(n • p)p 2 - 2(6 - 5i/)(n • p)(p • p) 
-2(1 - 2u)(n ■ r)(p • p) - 2(7 - 2v) [(n • p)(p • r) + (n • p)(p 
1 
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i(n • P)(n • r)p 2 (p • f) + 8(n • p) 3 (n • r)(p • r) 

+2(n • p) 2 p 2 r 2 + 5(p 2 ) 2 r 2 + (n • p) 4 r 2 - 5(n • r) 2 (p 2 ) 2 
-6(n • p) 2 (n • r) 2 p 2 - 5(n • p) 4 (n • r) 2 

-4(n.p) 2 (p.f) 2 -4p 2 (p.f) 2 ], 
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5(5 - 7z/)(n • p) 3 (n • r) + 10(3 - 5z/)(n • p) 2 (n • r) 2 

+3(17 - 35z^)(n • p)(n ■ r)p 2 - 28(3 - 8v)(n ■ p)(n • r)(p • 
-15(2 - 3i/)(n • p) 2 (p • f) + 2(24 - 77z/)(n • p) 2 r 2 
+2(9 - 29f)(n • r) 2 p 2 - 3(4 - 9^)p 2 (p • r) 

-2(12 - 37^)p 2 r 2 + 4(6 - 17v)(p ■ r) 2 j , 

D-3 = - — ^(3[927tt 2 - 10832 + 36 (48 - 5?r 2 ) v\ (n • p)(n • r) 
1536 L 

-6 [3(16 + 7T 2 ) - 2(45tt 2 - 464)^] (n • r) 2 
+ [11600 - 927tt 2 + 20(9tt 2 - m)v\ (p • r) 

+2[176 + 3tt 2 + 6(176 - 15tt 2 )i/] r 2 |. 



The last term in Eq. (0), SI (r, p), contains the parameters which parametrize the ambiguities of the 3PN Hamilto- 
nian. Its explicit form is 

2 

ft (r,p) = cj kinct ic[p 2 - 3(n • p) 2 ]^- + w s tatic^, (15) 

where akinetic parametrizes the kinetic ambiguity, and the static ambiguity is parametrized by w s tatic- 

The 3PN higher-order Hamiltonian H3PN (r,p,r,p), Eq. (||), can be reduced to a usual Hamiltonian depending 
only on r and p. Details of the reduction can be found in Sec. II of Ref. f|. 



III. THE AMBIGUOUS PART OF THE 3PN HAMILTONIAN 

It was found in Q|, [^) that the ambiguous part fl of the 3PN Hamiltonian, Eq. jl^), takes, in a non-center-of-mass 
reference frame and in non-reduced variables, the form (G, the Newtonian gravitational constant, is put equal to one) 

1 x ^ x f m a TO fc 

b 



fi(xi,X 2 ,Pl,p 2 ) = -jVV Static 



c 6 L^i \ static r A 
a b^a 



+ ^ kinctic ^^ [( Po ) 2 - 3(n ab • p a ) 2 ] } . (16) 

Z r ab ) 

In the following we shall show how center-of-mass considerations can impose restrictions onto the "kinetic" ambigu- 
ous term. To make contact with the intended center-of-mass motion calculations we transform ( |l6|) to the Lagrangean 
level. On this level, f2 appears as — fi with 



r!( Xl ,x 2 , Vl ,v 2 ) = 

a b^a ^ 



" l a , "'b 
^static 4 

r ab 



+ ^ ki netic ! ^ [K) 2 - 3(n afc v Q ) 2 ] } , (17) 
2 r ab J 

where the bodies' velocities are given by v a = p a /m a . 

In view of the relation between the total linear momentum P and the center-of-mass velocity V = X (X denotes 
the center-of-mass coordinate) demanded by global Lorentz invariance, 0, 

P = ^^' P = P 1+ P 2 ' ( 18 ) 

where H is the conserved total energy, H = Mc 2 + ff NR , the following expression (which is one of the contributions 
toP = dL/dv 1 + dL/dv 2 ) 

dfl dfl 1 \ - x m? m,}, 

D = — = -w kinctic — > > — n — v a - 3(n ab v a )n ah (19) 

a b=£a a " 



has either to be a total time derivative (in which case no restrictions are imposed on f2) or to be fixed by other terms 
in Eq. ([l8|) (note that no static ambiguity is involved here). 

The simplest way to show that D is not a total time derivative is to write it in the form 

n , u 2 . 2x T H , TOiTO 2 (mi-TO 2 ) - 1 

D = w kine ticmim 2 [(m 1 + m 2 )V J + — v J \d lj d li — , (20) 

M ri2 

(y l denotes the relative velocity, v l = v\ — u 2 , and here, MV l = m\v\ + to 2 d 2 holds) and to average D l over a circular 
orbit of the relative motion simply assuming that the center-of-mass motion is orthogonal to the relative motion. The 
result is (notice that the relative motion part in the Eq. (|2(]) is a total time derivative) 

< D >= -u kinctic 1 ^^(ml + mj)V. (21) 

r 12 

In the center-of-mass frame D % is always a total time derivative. This fits with Ref. Q where the kinetic ambiguity, in 
the center-of-mass frame, was shown to be influenced by coordinate transformations which additionally only influence, 
unimportant in this context, the static potential at 3PN order. 
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IV. CONCLUSIONS 



In this paper we have shown that the kinetic ambiguity obtained in binary point-mass calculations at the 3PN 
order should be getting fixed by center-of-mass motion considerations. The static ambiguity was known to be fixable 
by referring to e.g., the Brill-Lindquist (astatic = 0) or the Misner-Lindquist solution (ecstatic = — 1/8). The energy 
content of the latter solution is influenced by the topology of the involved non-simply connected 3-space. 
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